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A numerLcaL methocf for e-ciIvItis. Llcicar, two-cH mv.Ji b iona I clliptlL’ tniimdjjir'L' 
value profelecna ia presented. The method ia oa^enLlall^ tilt: Ritz procoilurt 

■JiLch uaej polytiornial spline fuiactioufl to appifoKLmate the enact eoluciour 
The spline functicms are constructed by defining ^ pol ynopilal function 
Over esoh of a sot of d isjoint EubduBneina. and Impo^ina 00rtain c^■^^^patlbll t ty 
c-onditiOnS along cOGimOn hoondaries between anbdoffl^ina . Tht main advantage- 
of ttie method iE chat it does not even require the continuity of the apllne 
fuhotiOna aefOsE chc boundaries between snbdcmaing. Tl-verefot# it ia eaay 
to COfiSCrott tlaSSOS of Eplino functions which will produce any specified 
rate of cortvergchcc. 


I* iNTRODyCTlOH 


J>LiriLi)^ tli& last fev years polyaonlal splines liave been u-^eJ to obt^ilri 

numerical soluticms of specific elliptic bcuiiddry value problems [ 2,3^43 

This paper prcseiita a general formulation of this type for solving two- 

dir^ensional elliptic boundary vaLv$ problomBr TtiO tfiftthod is osaentially 

the Rica procedu-ro applltd to a finite Jltufiicisloofil space of poLynoniial 

spLlfie functions. The formulation dan be Cutendccl to liijjhcr dilflensiono 1 
cas4&s and to slmultaneow$ systeas^ 

Ue denote by C S. bounded open domain In the Euclidean plane oich bouud.ary 
and cloeure G* Tor a function UeC^(G) ve denote by 1| 1*^111) G naximurn 

norm cmx max u(x) I and denote by Mu|[ the norm I ufioj^ds 

|(il=j e |o|^j 

Ol 

Here Q = (Ct^ , q:^) is the multiple index of letijjth .|df| » + CC^ and T is the 


eorrespondlng partial derivaCivOd We rlynote by |u ^ the scrai-nOttn 

^ 1/2 

f 2_/ iD^ul^ds') . Lot b = 12 b fx>D^ 

operator of oxdef Sin defined in {?. Suppose that b^c”(G) and G is c-f daaa 


be a real etrcngly elLiptio 


Let the bilinear form associated (jitb the operator L be given by 


2 


vhfirt + Then L - p'’^(h D^>. Wg aaauuic ttiat ?•(+>') 

|P 1^ 


satisfies Cllsj Carding's inequality tn (lie fallowing form: 


(1-2) pC;p*qp) [;£ ^WIm>G e ^ C-G') 

^ere G' la any oj^en Subset of G., 

For a given fgC^tG)j Itt UgC^tC) be the infinitely dlfferanc£gble solution 
satlefylng the e-quatlon Lu ■ £ and aaCiafylng sero D1 rich let d^jti on the 
boundary. Then 

<l.i> ^ 

Z. rnXJ-BIMEMSIOHAL SPLDIB llmcTlOWS 

Subdivide the domain G by d ftrid h into a finite number of suhdtuaalns 
Such that each aubdomain la a polygon infinitely differentiable curvilinear 

^ides+ We call any auth aubdofliain a ceU+ Let {S^.} be the set of gtid linea 
forming the sides of calls. We require that the Ifl.terlO'r of any gtid. line 
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muse nut uijntaln a veiftsK. sv-\i {Li-ngth of S^) • h, Wu denoto [C^} the? 

disjoint ijollcctlon of tntCir iufs o£ the cells and let G - ij G . 

U i t 

AssLime that ejeh grid line ia given by an equation of the form 


(L.4) ^^ = f. 


y - 3i (Si) 


lAere is the JistEitLce- along the grid Line 3^, measured from one o£ the 


vertices which it contiects. The unit tangent vector t^ ■ ^*'1 *^ 7 ^ ^ point 


on 


th^ grid line and the unit normal n^ “■ are defined by the equatlonj 


(1.5) 


U = 


df./d f. 

L. 


df. 


I 


t m 




r 1 

] 


2 lTi 


If u^^^ 6C™(G ), let denote the r"th derivative of in the 

1 n j r 
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dirtictitm cf fhe normal cif j defined if B ia one of LEie sidfs of G, . For ovory 

A. j 1 


cell, we define the following q^ianClClOs: 


Cl.6> 


“kr ^ ^k+1 I “nfr 


.CD 




defined for all It dfld c 
■whenever 

and O-^k-i-p-l-y- (m- 1 ) 


wfjSr-re p is some fixed integer 2 ;[) and indupepdent ofl-^. Also if w ia a vertex 


of a cell G. > dofitiC for 

1 Q Cc 


(wj for ell Of s'uc.Tl that fl:£| 0 f |5 m- 5 . 


iJe define a class of functions with dociaift G a* follows T Let 

deftot* the local coordinates of a cell C^* thtained Ijy tranelation 
of the origin to s point in G^. If then = u restricted to G. 

is given by an equation of the fom 


(l.D 


= a + a, x. + a , V. 
oo lo t fl-l ^1 


"d 

-. * + a ^ y . 
o,rt. i 


such that the folloHlng cotnpjitlbility oonditiona are satisfied: 


(i] For nil k and ^ and i^enever defined 
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«•« \i - ■ 

and = 0 le 3j » aCSfiaC^ 

(ii) At « vertex for gll |(S|£ rn-2 

(1.5') if V beLotifjs to both G. and G.± 

u L& 11 

and “0 if tho 

Vc iilao ttquite? that 
(1,10) rii + p 


Wow we impost the £oLlowing ermditianfl on the elass of eduniaglbLe: grids 
First we Lntrodyoe non-diriieflsiOfiitL toordinatea ■ y^/h,, 

and quantitlea ^ andU^^^= u^’^hl^l so that 


=■ A + A, X + 

DO lo i 


A y/ 

on^ t 


( 1.10 


+ 4 h 


6 


a.m vl.1 = J 3 ; 5 - {Ill'll 


- {hl'^ll3^>Li''(„) 


For B givpn cell let lie the vector space Epanned by the quantities 


,11 


* J '-fc 

and and let be the vector space spanned by the coafflcients 


Then ve hawe a linear CtanSf0rtnaEion 


(U13) Tj : 

We choose siifFlcltntlj' so Chat is PTl£o^ T^Ow let be the kernel 

of and let be the canplement of K^. Then T^ induces an 1 stmi-cirphiam 

(K14) t A. ~ V. 

' 111 

V'e reqult'S that i [a) VaTuee of are bounded unifonoly vith respect to all 
adtiiSiibIe grids atid the absolvts values of the eigenvalues Pf T. are'bouncledi 
flway from aero untfomiy with respect to ail admissible grids^ (b) There- CKiata 
a positive constant such that fPr any adinlsaible 


inf 

i 


length of S 

h ^ 


I 
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(e) There entsCa a positive cotiiitant timt ftJf cvf:ry arfmisEibLo grld^ 

and for e^ery ceLlj there exist# # &oirtt (k jV ) £n thiit c«ll sut-S that the 

y 0 

square {tiSntalned In the cell, (d) Modulo 

translfltiCJClf s rotation and scale factor, grid lines are re strict Ed to a 
finite numboi: of shapes, Ce) There erists a cons tan C such that the number 

of Cells thar any straight line through C tnccrsccts Is less than 

If the cells are triangular or reCtangulOP fa shape with atraight sides 

and if they satisfy conditions (b) and (c)* then it is possible to choose a 
Single integer n for n^ in Eq+ (1-^) such that condltlcn (a) Ig satisfied for 
gll auch celLa, To show this for triangular cells,, consider a cell in the 
form of sn tsoscelcS right-angled triangle with hypothenuse h units lonj and 
with origin of thei local coordinate syatem at one of the vertices. Since the 
cransformati-on for this particular cell doea not depend upon h, we can find 
an integer n for n^ in Eq. (L.7) such that conditicti {a) ig satisflE-d for all 
valuoa of h. Mow COnaider any other triangular cell which aatiafiea condi¬ 
tions (b) and (c). Eet n^, = n. Without Loss of generaTlty we can assume that 
haa tht origin ef the Local coordinate syattm at one of its vertices+ There 


s 


CJilstS a noiTi-singular lineau transformation S aicpping oi’ito 

specIficatIon of Che quantities and (l£ dcftncd) iniplles spQclficn- 

n ]r 

tion of all tiuantiCieE &£ thi: form 


1 

k-f-1 


f s’f 

J t,if’ 1 


O^r+r '■ 1 ^- 1 
f'>£> 

[j£.k^+ (iH- c' ) - (m-1) 


Where ±s I 
c j r 


n^^^ denotes the r'-tti derivative of In the direction of the 


tangents dofitved If la one of the aides of , 


Ther^:fore s induces an 1 aomorphlsm 


(1*17) 


^2 = t: v^, 


JiLso S Induce a an iaomorphism 


(i.LS) t Aj t; A^. 


Therefore wa have a -c ftranut a t i ve til a gram 


(1*19) 













1 , 


g 


Thlg , In Lurrii a ccTimutative tliajJtTaTn 


( 1 . 20 ) 




A., -+ V.. 


Therefore T. , is tionsingLiLar, Mow we claim that cPildiitiLHi; (a) ia aatlsfled. 

For, tf not, there exlats a sequence of triangular cells inch Chat che 

correaponling s&qutdtft o£ the LovesC abaolute eigenvalues of the tranaforr.e- 
Ctena T, tcncls to aero. EuC any triangular nell is ipetiffed by sIk quantities, 

jH 

namely the nOn'-Jimenatonal coordlnaCea nf ita veTtitC‘5+ Xhepefore the aequence 

of triangles can he cnbe^lded b* H tlosed enbapate Lying in the clostd annylu? 

in , lying between the unit sphere and the sphere of radius Jj- Since the 

annulus la cocjpactj thp; sequence has an accumulation paint. By continuity, 

the triangular cell G. corraaponding to the accumulation point aatlafies con- 

K o 


dlcltmg (b) and Ce> attd h*i singular T^^^ which contfadlcCe 






10 


In thfl sflM faahion, we emn th^Ose a ainele integer In Eq. U^ 7 ) whi.cli 
wtlJ Satisfy ccn-ditim (a) for all r&ttdngnlar cella provided that condtL>.flriE 
(b) and (c) arie satisfied. 

Tor practical applicationsj it Suffices tc aenmfcnC s -dfniieiiTi using Cri- 


angular and rectangular cells as well ns a finite number of specis1 ahapes to 
handle a curved houndary. 

The vector space "JtdJjfjp) is CDnstruct&d ^IS f^llovSi. Fflr «acVi cell C^s wC 
have an iaomorphism A.^ ~ Then a futictlon ue is defined by 

apecifying the quantities 


i 1 j i' j 

\t * \v 


^ = 4’" ' -a ■" 


for each grid line not 
in 


for each verteii not 
in 


and an eleoent of K 


Cl> 


for each eclL 


Lfe require that 

u^i. = OifSj isin^G 

and “'Oifw isin^G 


IL 


VJe BKtend the (semi) n-orms |' | ^ ] | ^ ] | > || ] - ] | ] bi 1 Linjii f f 

by restvicLing Che integrals to 

■5 "W 

LeC V be 5 ooetor ipacfi spimned by the quantitlea 
^G) with the baeie {e^}- Let E = @ K. wtth l«6i5 [J^-]h tct Ifl ^ 

L ^ 

we h^ve l^omorphisin 1 J ^(!f ^p) ~ Wr The liaagea of e. and 
a baala for iP) ■ Denote thia basla by ^“^^3 


L->™ B(' 1 -) to 


( 5 , and w not 

J 


V ® K. Then 


under I fore 


3 , FOEMlfLATIOM 


We seek a solotion such that the equation 


( 1 . 21 ) 




is aaCiafied for all pp)^ lie ahov Chat Eq. (1.^1) haa a unique solution 

and that Cn-^u) “ for 0'^,|a|!^-L . 

Tot a nuuietical solution of Eq. (1.2L) we reduce it to an equivalent matrix 
equation aa follows. We want to find v^W such that u ■ I ^(v) aatlaflea 
Eq^ (K21). Carrying out the integrations, “ (v^E v) where v = IC^j 


in 


u 


C d. ) t* thfi inner protJtiel; and E te a asatrix (B. .> where E * B(ii *u >■ Also 

ij ij 1 J 

f - f - 

where F^V isnd its i "th c tsmp ete nt la given by Thcrefoti^ 

■ (1.21) i$ equivalent to the cieErix equation.; 


(L.22) 


1 V ■ T 


If within each cell, wt appr-oxinat* the coefficients b and the function 

l-vp 

f by polynooiiala use of quadratufc formulaE fot integritlon of polynomial a over 
ElmplexsE f S] simplifies- the computation of B-. The computation is carriail -out 
by computing, Efu^tU^) tell by cell. 

4. COeii'EIlGENCE 

In the foil twins ^2' positive numberaT the choice of which 


does not depend uptn^ + 

Lemma 1: Let Uf'yfhjp)! tot be a gri-d line containedE in gGj.rQGj.j, and 
i i' 

let n • u reatrlcted to 0^ and u ■ n reatricted to . Then for |qi|^-1 




(1.23) 




u 


If Sj ■ then f&r |a|:^-^L 




sup 

(pi*y) e^j 




Proof; Suppoae that 3^ connects vC ft ices and WitSout lOsS o£ genavsHty 

WE? can assume that Is the origin and. that the K-aKls passes thmugh Let . 

X. be the length of S.i Waw consider a subapace the vector apace A 

} J ■ ® ^ 

spanned by '■■ % m-i^ ^ auhspace b-^ spennsd by 

I , * ^ . . {If m - 1, then (v ia one- 

t OS J aUm-2 t a ^|cf]^-£ Ojm-L o 

dimenslonaL, spanned by Let j : Q A be the inclusion map Add 

p : -f ( 5 ,^ he the projection mapr then It Is easy to checlE that the map 

by the tCnnpoS 1 Cic-n poTj,oj: la an isOniecphism. and therefore 
CA^. Moreover aincO jrld Lines are reaCricCed to a finite nuniher of 
shapes and /h ^ tt followa that the absolute VAluee of tJr- eigon- 

va^ues of are bounded away from zero uniformly with xepect to -Jjf - 

Let be the cOmpItment of in and 01^ be the complement of 
in (^.^p [)ij be the corre-spondlng spaces defined for the cell 
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Lut aciif Ipt anJ VifK he tin' vi'i-ti>r3 c pt re sptind i Jiji ti> i\. 

Let the pTC:jetl:lon of u onto be v iintl l&t v = jj.^ 

and Vj^E . Let vnique po tj-ndenf a I c orre appmiit^ & to 

beP=e +E^?: + **,+G ,’^'1^'^. SimtlarLy^ tbe corresponding v^sCCor 

pD ^ 10 ± o 5 m- L 1 

v^e detetminffi ft unique pc-l^niMiiial ?* = “‘‘^o ,[n-l'’^i ' 


V = . Therefcr* P = P * 

o 0 


Let d = T^^[<ocBVj^)) and d - tJberC ■ 


__ *■ * 

Let reatricted to L^.^* q = ^±1^ Then poT^ + poT^ is the aero ciap and 


Tb ‘i --q 


poT* = T^* Therefore Ct^ Cd^>) ^^Ich ImplleE that if d^ - 


|Ct|<JH. 


nr 


We define the foLLowing nonas: 

?or q - q^^, ..^. the norm ||q||^^^ ^ 


Fciif r 



^lo' 


,t )c'Cti define the norm 


UM !*■!'=’', 

(1 hj J 



For S = Cs 

]TI p O 


s define the norefl ||e||= ma?; ]*j^j | 

’ i (i?j} 
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For dcfina tine noriu '| |s [ [ " \\^\ \ I 1 t' 

Tlieii 1 jdp] laJCj l 1^*1 I I ■ Similarly [ [^*^3 | I^J M ^ 

T.[((a + d ) $ projection of k onto K^ the vector 

k - (k , ... k > and let d £ k * f. and 6l ® k' = {'‘. Then: 

^ 1 Ot* Ojtlj i 1 L J. i * 










Therefore aup [D^ ^(K,y)J^-i|r ( j | |j. | | + j ^ 

(a:.y)eSj 111 I 


Moi* we claim that 




([{.IP I “I 


To proTO this obfi-erve thiit 


, m-1 
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The quaJeaCic t ■ ] 


■[E 


Inf 

0^2 g 


ia poBitlve definite an the liti&ar space ^ Thc-tefore 

1/2 

— ^ [ ^ 1 ^ 1 = X ^0. \ is s funeti'£>n of 1 and h ■ Suppose Inf ^ * 0. 

]]i|] 1 i- ^ ^ 

Then there eKlBta a sequence or cells such that the c or re Sp-an-d Ing 

aeqnetlCt 0 * Non each Cj^ contains a paint the aqusre ■ ■ 

-at - ^ the blllneat fans [ , ] 

xestrlctEd to the square Then there exists a sequBiwe of point* * 

m 

l\ ® earreaponding sequence 5 uCti that and 

■s’ . 

0. By cocrpactness, there eKlste a subsequenac ^ 

K n J rJ J 

verging to (X^ ) U fiath that ||z^J| - 1 and ~ ^ 

■0^4* 

which is -A contract let Ion of the poaitlve-deftnitatiess of [ , I’. This proves 

the clain* Inequality U-^i> follows. The proof of Inequality (l.Z4> 1* slmlUr* 
LemiM 3: Let Then 


(1.2^) 
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Proof! 


Hilhout Lobs of *eii€rallty we nay iisstHiiE r.bat none 


>f tM lifLOi 


Is para 


lUl to tho coordinate ^ funCtioo o’ 0" 


Tif 


by loCCioe 


u' = u on G atid u‘ = 0 elaevhere* For flny Une ^ Lnter 0 cCClTi& G^. P 


pcinta 


.yj C-b-V,) i" 3^ air points h.vlo£ coordlr..os 


(ic,y^) irid lying In G be between and ^ . 

I ■: . ■■ -’' 

"Tfl-l. deilne a. fonctlon 

J ^(0 -(!>„«'<-,*' C*b-.> ■’'c> ■='<.> =«■. 

0 othetiflse 


Ttien 


J (o t= pUc,»Ub wntibbbu, »nd bv.r op™ «gn««* 


ftefitie ° 

and 




dj 


a 


Let J^tt) 


t ^ t^ 

C-E:j <0 

f„ ( in <mb b£ tl» bfs" s6S».™t» listed dbc«. thbn 



K+l 


fi>r ,■■ Fcr any pfllPt 


ta{tj)p 


JJO - Jj^(t 


i+)+ 


(t) 




t 

‘=[2i ‘VV>-VV>’]" ^ = [1 

^ j-" ■ "o' 


* 2i[|; fV^) ^ "'[/ '•'» *’ 


)|^ dT 


by CiU-Cjiy-'Schwafa Intqwality + From and. (1.2^) , 


2^£j_^Ct^+> - Jg(tj-)J^ ^ ^ ^Sl“L,e^ 

Fv^ Inequality 


1 <G, 


If 




J <. 
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i. . • 




Thertfo 


rE 


<!’ * hj |u|f „ /li^ < ^ |ii|^ 


k "’°k 


h ' 


Condi cion (a)' an the claaa of adiElsalble grids Intpllea that i s- ^ 

h 


and 


compactneaB of Implies that t f 


Therefore 


h ni,G^ 


«n^ |||D^u(||^^ ^ ^ _§ M 

fly aucceaaive integrations In a sltailar fashion, 


for [ml ■ Di-1 


I ^O.G^ 


'3 n 


m 


for o S; n-1. 


The lenime folLowik 


Va state vlthouC proof tho follovlng special forn of Che Green's 


formula. The ptPOf is slfliiLar tO the one given In reference w 


V- 


Lemiiia 3. There linc^ir differential opcretcirs tor Ck, y>^ 1^50^ 

„Uli bonded unifL^efflly with fL^spect to isll U,y) =3ml allk, su.h 

that S', is of ordei: and for all ueC '(C]l and ' 


( 1 . 2 ^) 


f jjlju “ ^tm ill) "I" 5^ r ^ ^j] i 

^ ^ S^4f.G j-0 \ 




n,J 


Kl”! 







whtrt the grid lin« S, = aG.OSG,' or SCH BG. 

^ k 


Lciiiiaft 

For a jiven u e: 

C“'(G> 



(.1. ^ Ui - 

G 


for 

all ^ e^tKp^ iJinJ all 1^ 

Pronft 

If a = (aj^» 

^ lot Cd 


1.1^ 1 ari(i 


U, yf 

= X y 

Ihcfl hy Taylprs' cheoreTH 


" E . 

\a [ fp+m 




OL «2 ^ ^ 

(y “ ■*' ft^Kjy) 


uhere 


2 L 


llllf RtK,y,>||| 

0 


and 


d = diet an cc between polTits tTt,y) and y^>. 

liDW any grid, Line Sj^ is given by the eqwatipfis ?£ * f|^ ^ ^ 

•30 

Let y^} e point in S^. Since Sfid are C functions and since 

the nunib-dr o£ shapes of grid lines is finite. It Is possibls to choose 
snoh that 


(1.29) 


!■ .n 
n* j 


p+m-J 

r-^0 




w 


where I I I I I I ^ 


p+ra-j+l 

^10 


Substituting 


(1*^9) In (1.27) 


35 


f 


- 3( 0,u} I ^ 


sra “'r • 




j=0 r*CJ 


f L i * 

{ (D ,<» . + Ah. 

J a,i nj ■ J I 


m-l p+J.*l-it . 1 ^ 

E [ £{ E 


S^eiG J-V fO 


+ 


/ 'V}] 


"l>*» I I Ir^^j {«^) 1 i I « C,(j 6 f 


Sirica I aatleflas ch« Compatibility C0ll41 tlutis 


1, b i:b 

^1 rj 


if 


B - 9G nw ( 
\ \ 


■y 

=” 


if 


ac. nsG 


&y Ceuchy^SchVATz Inequality^ 


23 







[,SJ. 




{■u&lna InequaLtptes (1.13) ^ 


^ 'lO 




m ,C_. 


QJ L 


^l(p+j+2-Tn>j^ 


(ajffia o£ G} - 


f Si 




Slraiiaf iy 




B. eSC s, 
k H 




2 ( p + l ),,|2 


iii>G 




Therefore 
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1; Eqiiaticm £1+2:1) has a unique solution Ug >C^jp) atid 


ti.il) 


-u 




p+^ 


(1.31) 


1u-ui ^ < C h 

^^G^ ^ ^27 


P+1 


im I 


FfOOft Suppose for ^ ECfpt^S-) “ ^ tp' 


** I I |“l I L-1 G “= <) tiy Inequality (1.2^) 

+ ■. u = 0 

i' 

Eq. (l.ll) has a unique Ue?^(^ip). 

Using Taylor's theoreai;^ within each cell* we e@U write the exact boIu- 

*1 '^1^ ^ polynoiiLia 1 of degree 

*1 for 


- U^>“ 
V a 


u,+l 


tioci 

as U * 

'u(K^,yj^) + R(iCj.*y^) whera 


rt 

V ^ 

Let 

’■<1 = 

U^J-'ujJ .»d 

iriT^ 

a 

Then 

i""ii 

1 * 

l,.^,i.., 


23 


Lee <"Lf 


I ( all 




_-l 

Let a" T, (d^“(A i • ^ * A ^ 3 + 
1 oo &ttl. 

£ 


Then ]A^j| 


„ . fl £+1 

‘^U '' 


Ut "u(x^.y^) • + dj, + ... 


"1 


v+ier£ a, , = A /h^'''J 

1 > j * J 


L-ki- - 1g| 


IP„"'’IL,o^’=“1? " ‘ 


for a .^n 


Let '^■u "■ 'u + "u and u ■ u - u , T^en 


and 

M ”g~ II 

1 f, -SC., , 

1 P i L-H 

]“i 

Let 

ii 

^ u - 



f- B 

ECa^ju) 

p 

(Li) 



= B 4 { 

/CAh) {LuJ - 

By Tne< 5 T.La Lit y 

a.ag-) 



- ^ CiLi^u} 


p4l 
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Nov the Cs-ueh^-SchwdTS InequaliCy implies thet 




cp+l> 

Therefore j'BC^u^u) h c InequalltieQ (l.tOJ and (1. 26). 

Inequalities <1-3) and (1,33) imply that 

(p+1) 

: l^^lu*C ^24^ 

? ■.. ■ O 

’■ S 

Inequality (1.26) implies that 




'■-IIL.i.Vl!l^“-IIU-i,o/lll^"IIU.«^ * IIMIU., 


i:- . O 


tp+U 
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